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What information can we obtain from the 
yield ratio n~ /n + in heavy-ion collisions ? 
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Abstract 

The recently reported data on the yield ratio ir~ /ir + in central 
rapidity region of heavy-ion collisions are analyzed by theoretical for- 
\Q [ mula which accounts for Coulomb interaction between central charged 

fragment (CCF) consisting of nearly stopped nucleons with effective 
charge Z e g and charged pions produced in the same region of the 



phase space. The Coulomb wave function method is used instead of 
the usual Gamow factor in order to account for the finite production 
range of pions, j3. For Gaussian shape of the pion production sources 
it results in a quasi-scaling in (3 and Z e g which makes determination 
of parameters (3 and Z c $ from the existing experimental data difficult. 
Only sufficiently accurate data taken in the extreme small m T -m n re- 



gion, where this quasi-scaling is broken, could be used for this purpose. 
PACS numbers: 25.70.Np, 03.65.Ca 



The ratios of pionic yields tt~/tt + in the central rapidity region in high 
energy heavy-ion collisions have been recently reported by E866 |lj and by 
NA44 p| Collaborations. In both cases a strong excess of yield of n~ over 
that of 7r + in the low transverse mass region has been found. Although 
its theoretical origin is still point of debate, the most obvious reason to be 
checked at first instance is the effect of the Coulomb final state interaction 
taking place between charged pions and central charged fragment (CCF) 
consisting of nearly stopped nucleons with effective charge Z e s and produced 
in the central rapidity region ||, cf. Fig. 1. In the approximation neglecting 



the possible effect of the finite range of the production of pions as seen, for 
example, in Bose-Einstein correlation experiments, the observed results are 
given by the ratio of Gamow factors only : 

^ + ( Prcl )/iV^(p rcl ) = G( V )/G(-rj), (1) 

where Gamow factor G(rj) is defined as G{rj) = 2nr] /(exp(27r?7) — 1) and 
1] = Z e g m n a/ p rel . Here a, and p re i are the fine structure constant, mass 
of 7r meson and its relative momentum in the two-body system (7T-CCF), 
respectively. The energies considered here are high enough to neglect the 
influences of charged spectators produced in the fragmentation regions, cf. 
Fig. 1. 

Such problems have been already considered before but either at much 
lower energies |4], [| or with the use of Gamow factor only || . In this work we 
shall analyse high energy data and use the Coulomb wave function method 
as developed in instead of simple Gamow factors to describe the final 
state Coulomb interactions. This allows us to account for the finite size of 
the emitting source (or, equivalently, for the finiteness of the pion produc- 
tion range) characterized by parameter f3. The question we are interested in 
is: can one obtain from the observed yields of n~ /n + production mentioned 
before any valuable information about both the effective charge of CCF, Z e g, 
and the range of interaction parameter (3 ? We derive theoretical formula 
for the pion production yield based on the Coulomb wave function convo- 
luted with some pionic source function p(r) and apply it to the analysis of 
recent experimental data on Au+Au collisions at 11.0 GeV/nucleon M and 
on Pb+Pb Collisions at 158 GeV/nucleon 0. The emerging valley-like struc- 
tures in the map of the x 2_va hies made in f3-Z e ^ parameter space is shown 
to be connected with the quasi-scaling behaviour of integrals of the square 
of Coulomb wave function convoluted with pion production source function 
||. This finding shows that the expected determination of parameters (3 and 
Z e Q from the existing experimental data is difficult if not impossible. 

Let us consider Coulomb interaction between 7r + with lab rapidity and 
CCF of some effective mass M e g and with rapidity Yqm (which we assume 
to be equal to the cms rapidity of the colliding system). In the rest frame of 
CCF the 7r + momentum p„- is given by 

Pvr = ( P L , = ( m-rsinh^ - 1cm) , Pt ) (2) 

where m T = \jm\ + p| and p T is its transverse momentum. In this frame 
the Schrodinger equation for the relative motion of the two-body system 



(7r + -CCF) is given by 



2 7 p i 



Pr ^eff 

2/i r 



V» r (r) = E r ip r (r) 



(3) 



where r is the relative coordinate of tt + in the rest frame of the CCF, 
p r = —ihV r and p is the reduced mass of our two-body system: \i = 
m w M e ff/(M e g + m n ) ~ 0. One can find that the wave function ip r 
is given by the following confluent hypergeometric function ]7|, |10 1 : 



^r(Pr,r) = 

T(l + tr])e- nv/2 e iPr - r F( - irj, 1, i (p r r - p r • r)) 



(4) 



where p r = M e s p 7r /(M e fj + m n ). Convoluting now eq.([|) with some source 
function p(r) (i.e., distorting the Coulomb wave function accordingly in order 
to account for the finite production range of pions), we can now calculate 
single particle spectra of produced 7r + 's: 

N* + (p r ; V+ ,f3)= f dhp(r) Vr(Pr,r 



oo oo 



xI R (n,m) (2 p r )™ + "\ 



(5) 



where A n {jj + ) = T(ir] + + n) / {T{irj + ){n\) 2 ) and r/ + = +Z e g pa/p r . Assuming 
now that p{r) is given by Gaussian distribution: p(r) = ( ^jjgjf ) ex P (^) > 
in which case 



/"OO 

I R {n,m) = 4tt / drr n+m+2 p(r) 
Jo 



2 ^ +mr /n + m + 3 



(cf. Aic in ref.0), and using the following decomposition of Q 2 : 

Q 2 = p 2 L + pi = p 2 L + (m T - m^) 2 + 2m 7r (m T - m % ), 
we finally obtain that 

N w+ (m t - rj + ,P) 



oo oo 



fixed ra + m + 1 



v(n+m)/2 



x (Pl + ( m T - ^tt) 2 + 2m 7r (m T - m^) 



(6) 



For the ir~ production case the corresponding N w yield can be obtained 
by simply changing the sign of the Sommerfeld parameter in eq.@: r] + — > 
—7] + = T]_ = —Z e g fia/Q. Therefore theoretical formula for the ratio of the 
production yields we are looking for is given by: 



0. For Au+Au collisions data Ycm=1-58 and y 7T =2.10 in eq.(0) whereas for 
data on Pb+Pb collisions Y"cm=2.90 and for y n we use the averaged value 
over rapidities of 7r employed in the respective data analysis. Applying the 



one discovers (after using a lot of CPU-time) the valley-like structures in the 
maps of the x 2 -values made in (3-Z e g parameter space, cf. Figs. 2a and 2b for 
E866 and NA44 data, respectively |TT| . The values of \ 2 along these valleys 
are almost constants and equal to x 2 ~ 35/30 and x 2 ~ 52/53 for E866 and 
NA44 data, respectively. The traces of the minimum x 2 -values are almost 
linear. 

These results strongly suggest the following quasi-scaling behaviour being 
present in eq.(|]): 



where A>0. Fig. 3 showing the results of calculations of eq.(|8]) fully confirms 
this supposition: both for N n+ and N n and also for their ratio N 71 /A^ 7r+ . It 
means, and it is demonstrated in Fig. 4 for both sets of data analysed here, 
that present experimental data on yields n" /it + are not able to determine pa- 
rameters (j3, Z eS ) uniquely. As one can see we can explain E866 data by two 
sets of parameters: (/3, Z e g) = (1.0 fm, 24) and (3.0 fm,72). Similarly NA44 
data can be also described by two sets of parameters: (/?, Z e ff) = (2.5 fm, 40) 
and (5.0 fm,80). Only at extreme small m T — m n region, where such quasi- 
scaling in eq.(|8]) seems to be violated (cf. Fig. 5), such unique determination 
could be (in principle) possibly achieved (but only for sufficiently accurate 
data). 




minimum x 



2 -fitting method when comparing eq.(|7]) with experimental data 



N w (m T — m^; XxZ cS , Ax/5) 
w N* + {m T -m^ Z eS , (3), 



(8) 



Summarizing: we have analyzed recent high energy data for yield ra- 
tios 7r~/7r + measured in the central rapidity regions using our theoretical 



formula ([F[) calculated by convoluting the square of Coulomb wave function 
with Gaussian source function for pion production. In this way we have ac- 
counted for the distortion of Coulomb wave function caused by the finite size 
of the production region. With this formula we find the valley-like struc- 
tures for the maps of x 2 -values in the /3-Z e ff parameter space which can be 
attributed to the quasi-scaling behaviour of eq.(H) @- Results obtained by 
using only Gamow factor || do not show this property. It means, as was 
show in Fig. 4, that it is difficult, if not impossible, to determine parameters 
(3 and Z c s uniquely in this case. They could be determined, if at all, only 
by using very accurate data obtained in the extreme small m T — region 
where this quasi-scaling property in eq.(|8p is broken (cf. Fig. 5). 
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Figure Captions. 

Fig. 1. The picture of collision: charged 7r's and CCF are produced in the 
central rapidity region whereas spectators populate both target and 
projectile fragmentation regions. The transverse energy spectra of the 
charged 7r's are affected by their Coulomb interaction with the CCF. 

Fig. 2. (a) The valley-like structure exhibited by values of x 2 ( om y numer- 
ators are shown, denominators, denoting the NDF's, are equal to 30 
here) in the j3-Z e s parameters space for E866 data [|J fitted by using 
eq.(0). The minimum x 2 - value is about 35/30. (6) The same as for (a) 
but for NA44 data @ (here the NDF's are equal to 53). The minimum 
X 2 - value is about 52/53. In both cases the Z e s from fits using Gamow 
factor only are shown by the star (Z e g =7.3 for x 2 — 44/31 for (a) 
and Z e ff = 6.9 for x 2 — 164/54 for (b)). For E866 data eq.(^) has been 
renormalized accordingly ||11|| . 

Fig. 3. The quasi-scaling property of eq.(|6|) for: (a) N n and N n+ pion 
production yields and (b) for their ratio N n+ /N n . Three different 
(scaled by factors A = 2 and 3) sets of parameters: (/?, Z e s) =(1.0 fm, 
24), (2.0 fm, 48) and (3.0 fm, 72) give approximately the same results 
(solid line). Dashed curves represent the results obtained using Gamow 
factor. 



Fig. 4. (a) Comparison of eq.(|7|) with E866 data [[|] using parameter sets 
(/?, Z cS ) = (1.0 fm, 24) and (3.0 fm, 72). (6) The same for NA44 data 
using parameter sets (/?, Z eff ) = (2.5 fm, 40) and (5.0 fm, 80). For E866 
data eq.([7D has been renormalized accordingly ||11|| . 

Fig. 5. Example of violation of the quasi-scaling property of eq.@ for small 
values of the variable (m T — m n ) shown for N w+ /N w and compared to 
NA44 data § (here (p, Z e{i ) = (1.25fm,20), (2.5fm,40) and (5.0fm,80)). 
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